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INTRODUCTION 


Ocular  trauma  during  military  conflicts  has  steadily  increased  from  0.5% 
in  the  civil  war  to  13%  in  present  day.  This  increase  is  likely  associated  with 
the  advancement  of  weaponry  and  the  increased  use  of  explosive  devices. 
The  majority  of  eye  injuries  from  an  explosion  can  be  classified  as  either 
open  globe  or  closed  globe.  Open  globe  injury  is  often  readily  identifiable  and 
typically  undergoes  urgent  surgical  repair.  However,  closed  globe  injury  may 
not  be  detected  immediately  and  can  result  in  a  series  of  sequelae  that  lead 
to  visual  dysfunction  months  after  the  blast.  The  progression  of  closed  globe 
eye  injury  and  visual  degradation  following  blast  exposure  has  not  been  well 
characterized.  Furthermore,  it  is  unknown  if  there  are  early  indicators  that 
denote  an  increased  risk  for  developing  visual  dysfunction  following  blast 
exposure.  Therefore,  the  objectives  of  this  proposal  are  to  investigate  the 
temporal  progression  of  eye  injury  from  blast  exposure  and  identify  early 
predictors  of  visual  dysfunction.  We  propose  to  accomplish  these  objectives 
by  first  identifying  the  probability  of  military  personnel  developing  visual 
system  injury  after  blast  exposure,  and  determining  the  time  point  after 
blast  exposure  that  visual  system  injury  becomes  identifiable.  Next,  we 
propose  to  systematically  evaluate  the  time  course  of  visual  system  injury 
from  blast  exposure  using  our  existing  rat  model  for  blast  traumatic  brain 
injury.  From  these  experimental  studies  we  can  identify  early  predictors  of 
visual  dysfunction.  Finally,  we  will  evaluate  these  early  predictors  in  a  clinic 
setting  to  verify  their  usefulness  in  real-world  scenarios.  By  understanding 
the  temporal  and  chemical  progression  of  eye  injury  from  blast  exposure,  we 
can  establish  early  identifiers  of  visual  system  injury.  This  will  enhance  our 
diagnostic  capabilities  and  lead  to  the  development  of  time-dependent 
treatment  strategies  to  mitigate  the  loss  of  vision  in  military  personnel. 
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BODY 


Aim  1:  Investigate  the  progression  of  visual  system  injury  in  service 
members  exposed  to  a  blast. 

SOU/  1:  Retrospective  analysis  of  the  military  personnel  with  visual  system 
injury  attributed  to  blast  exposure  but  not  immediately  identifiable  at  the 
time  of  the  blast 

Obtaining  IRB  approval  from  all  parties  was  an  unexpectedly  long 
process.  The  University  of  Utah  IRB  approval  was  obtained  Oct.  12,  2012. 
The  VA  Hospital  of  Salt  Lake  City  IRB  approval  was  obtained  Mar  28  2013. 
The  HRPO  IRB  approval  was  obtained  June  28,  2013.  To  speed  up  the 
subsequent  tasks  for  this  aim,  a  project  coordinator  and  medical  student 
were  hired.  Since  all  approvals  were  obtained,  a  formal  request  for  MR#s 
has  been  made  to  the  University  of  Utah  Hospital.  Results  from  the  query 
are  currently  being  evaluated.  Accessing  veteran  hospital  records  proved  to 
be  more  challenging.  To  overcome  this  issue,  the  PI  recently  met  with  the 
polytrauma  group  at  the  VASLC.  A  collaboration  was  established  to  identify 
blast  victims  from  their  databases.  Current  security  approvals  are  being 
sought  and  record  evaluation  will  begin  late  October.  It  is  anticipated  that 
the  study  will  only  be  delayed  by  3-5  months. 

Aim  2:  Investigate  the  progression  of  visual  system  injury  following 
blast  exposure  in  an  animal  model  and  identify  early  indicators  of 
visual  dysfunction. 

SO  1/1/  1:  Manufacture  new  blast  device  with  enhanced  mach  stem 
overpressure  capabilities. 

The  shock  tube  for  the  experimental  studies  was  purchased  from 
AgileNano  (Fig  la-b).  It  has  a  6"  internal  diameter  and  15'  length,  which  is 
sufficient  for  blast  testing  in  the  rat  model.  The  shock  tube  is  instrumented 
with  side-on  mounted  PCB  113B26  pressure  transducers  (PCB  Piezotronics, 
NY)  installed  every  two  feet  along  the  driven  section  (Fig  Ic).  These  sensors 
can  measure  up  to  500  psi  and  are  monitored  at  1  MS/s  via  Labview 
(National  Instruments,  TX).  Extensive  work  was  done  to  set  up  and  validate 
the  behavior  of  the  shock  tube.  This  included  characterizing  the  shock  tube 
across  a  range  of  driver  pressures  and  membrane  thicknesses  (Fig. 2). 
Additionally,  work  was  done  to  devise  an  appropriate  filtering  system  to 
acquire  clean  pressure  data  from  the  transducers.  After  extensive  filter 
design  using  spectral  analysis,  a  cutoff  frequency  of  275  kHz  was  selected. 
Upon  data  collection,  the  pressure  signal  passed  through  an  analog  anti¬ 
aliasing  low  pass  filter  (Fig.  3).  A  digital  50  kHz  filter  was  then  implemented 
post-hoc  in  MATLAB.  This  protocol  was  found  to  eliminate  the  majority  of  the 
noise  without  significantly  decreasing  the  peaks  of  the  signal. 
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Figure  1.  (A)  Shock  tube  purchased  for  blast  studies.  (B)  Close-up  of  the  driver  section  of 
the  shock  tube.  (C)  Overpressure  recorded  at  seven  locations  along  driven  section 


Figure  2.  A  variety  of  membrane  thicknesses  were  evaluated  to  determine  the  appropriate 
experimental  conditions  for  the  desired  peak  overpressures. 
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Characterization  also  included  creation  of  a  finite  element  model  of  the 
shock  tube  to  examine  the  behavior  of  the  shock  tube  with  different  driver 
gases  and  to  determine  the  appropriate  location  for  rate  placement 
(Appendix  A).  We  also  verified  that  the  blast  wave  produced  by  the  shock 
tube  appropriately  modeled  the  Friedlander  wave,  which  is  characteristic  of 
open-field  blasts. 

A  decision  was  made  to  place  the  animal  inside  of  the  shock  tube  rather 
than  outside  of  the  muzzle.  This  allows  the  blast  exposure  to  have  a  duration 
of  approximately  2-8  msec.  It  also  reduces  the  variability  in  shock  tube 
strength  inherent  to  placement  outside  the  shock  tube.  A  rat  holder  was 
designed  and  3-D  printed  to  allow  for  a  test  animal  to  be  placed  securely 
within  the  blast  tube  and  subjected  to  a  side-on  blast  (Fig.  4).  The  primary 
design  consideration  was  allowing  blast  exposure  to  the  head  while 
preventing  blast  trauma  to  the  torso  as  lung  damage  may  be  fatal  at  the 
blast  levels  currently  being  used.  The  protective  measures  include  placing 
the  body  of  the  animal  outside  the  path  of  the  blast  wave  and  wrapping  the 
body  in  a  Kevlar  vest  for  further  protection.  Initial  animal  experiments  have 
not  shown  any  signs  of  lung  damage  so  these  measures  have  been  deemed 
successful.  Low  vibration  of  the  mount  was  verified  by  placing  linear 
accelerometers  onto  the  mount  and  measuring  motion  during  a  blast  with  no 
animal  in  the  tube.  A  maximum  displacement  of  <0.5  mm  was  recorded. 

To  ensure  the  safety  of  the  operators  of  the  device,  ear  protection  is  worn 
during  use  of  the  shock  tube.  However,  a  recent  inspection  found  that 
decibel  levels  outside  of  the  blast  room  were  too  high.  We  are  currently  in 
the  process  of  manufacturing  a  sound  absorption  chamber  for  the  end  of  the 
shock  tube  (Fig.  5).  Testing  is  underway  to  ensure  sound  outside  of  the  blast 
room  is  reduced  and  the  pressure  profile  is  unaffected. 
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Figure  4.  (A)  Animal  mount  designed  to  place  the  head  inside  the  shock  tube  and  position  it 
for  a  lateral  blast  exposure.  (B)  The  end  of  the  mount  is  clear  in  order  to  monitor  the 
respiration  of  the  rat  during  the  blast. 


Figure  5.  A  sound  absorption  chamber  was  designed  and  installed  to  the  end  of  the  shock 
tube  (A)  to  decrease  the  decibel  level  of  the  shock  tube.  (B)  The  chamber  is  lined  with  sound 
absorbing  foam. 

SOl/1/  2;  Establish  randomized  testing  protocol  for  12  experimental  groups:  4 
severities  of  blast  exposure  (140kPa,  260kPa,  450kPa,  and  control)  each 
with  3  survival  periods  (1  day,  1  week,  4  weeks). 

Post-injury  time  following  blast  exposure  is  the  most  important  factor  of 
the  study  and  new  literature  suggested  that  4  weeks  may  not  be  long 
enough  to  explore.  Therefore,  we  increased  the  number  of  survival  periods 
to  4  (1  day,  2  weeks,  4  weeks,  8  weeks)  and  decreased  the  number  of  blast 
levels  explored  to  2  (control,  225  kPa  and  450  kPa).  We  have  already  begun 
to  see  some  ocular  injury  at  the  low  blast  level. 

SOl/1/  3:  Complete  experimental  studies  investigating  temporal  changes  of 
visual  system  injury  from  blast  exposure.  (Year  1-Year  3) 

To  verify  that  our  experimental  design  (Fig.  6)  will  result  in  a  low 
mortality  rate,  the  low  blast  level  is  being  evaluated  first  before  proceeding 
to  the  high  blast  level.  Briefly,  animals  are  administered  carprofen  one  day 
before  the  blast  for  pain  management.  A  baseline  of  vision  functionality  is 
established  before  the  blast  using  the  metrics  detailed  in  SOW  4.  On  the  day 
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of  the  blast,  the  animal  is  anesthetized  using  ketamine/xylazine  or 
ketamine/dexmetomedine  and  placed  in  the  holder  within  the  shock  tube.  A 
pressure  wave  is  initiated.  The  animal  is  removed  from  the  device,  allowed 
to  recover  from  the  anesthesia,  and  then  returned  to  the  animal  facility. 
While  animals  do  not  show  signs  of  pain  following  the  blast  exposure, 
carprofen  is  administered  the  next  day  as  a  precaution.  The  vision  metrics 
are  then  repeated  the  day  after  the  blast  and  every  subsequent  week 
following  the  blast  until  sacrifice  (see  SOW  4).  At  sacrifice,  the  eyes  and 
brain  are  harvested  for  later  analysis  (see  SOW  5).  To  date,  5  animals  have 
been  evaluated  (1  2-week  survival,  2  8-week  survivals,  and  2  blast 
characterization).  It  is  anticipated  that  40  animals  will  be  evaluated  by  Dec 
2013. 
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Figure  6.  Procedural  timeline  for  blast  animal  studies  and  vision  testing. 


Two  animals  were  used  to  determine  the  effect  of  the  blast  on  intraocular 
pressure  (lOP).  The  lOP  sensors  used  for  these  studies  were  damaged  due 
to  the  high  pressure  loads.  We  are  currently  reviewing  new  pressure  systems 
to  continue  these  characterizations. 


SOW  4:  Complete  animal  OCT  and  behavior  studies  to  test  visual  function  on 
all  animals  prior  to  blast  exposure  and  at  1  day,  1  week,  or  4  weeks  post 
injury.  (Year  1-Year  3) 

Behavior  Testing:  Optokinetic  Tracking  Test 

In  Year  1,  the  behavioral  test  system  was  designed,  built,  and 
programmed.  The  system  is  made  up  of  four  television  monitors  to  display 
the  visual  stimulus  to  the  animal  in  360  degrees  of  rotation  (Fig.  7a) 
Behavior  is  monitored  using  a  closed  circuit  video  system.  The  coding  to 
control  and  display  the  behavioral  stimulus  was  created  in-house  using 
MATLAB  software  (Mathworks,  MA)  with  the  psychtoolbox  software  package. 
The  output  of  the  behavioral  system  is  an  estimate  of  the  visual  acuity 
threshold  from  zero  (best)  to  one  (worst).  The  relative  visual  acuity  is 
tracked  week  to  week  to  identify  visual  degradation  (Fig.  7b).  For  increased 
accuracy,  each  animal  is  tested  three  times  at  each  testing  day  so  that  an 
averaged  acuity  may  be  calculated  with  measures  of  uncertainty. 
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Figure  7.  An  optokinetic  system  was  buiit  to  test  the  visuai  change  in  the  rats  over  time.  (A) 
The  system  consists  of  4  LCD  screens  which  vary  the  biack  and  white  contrast  bars 
throughout  the  test.  (B)  Visuai  acuity  is  assessed  by  the  rat's  abiiity  to  see  the  distinction 
between  the  biack  and  white  bars.  The  program  automaticaiiy  adjusts  the  contrast  untii  an 
equiiibrium  is  reached. 


Optical  Coherence  Tomography  (OCTj 

OCT  can  be  used  to  image  the  retinal  and  evaluate  changes  in  retinal 
thickness.  A  thinning  of  this  layer  is  often  reported  to  be  associated  with 
visual  dysfunction.  Furthermore,  histological  findings  such  as  retinal 
hemorrhages  or  retinal  detachments  can  be  detected  using  this  technique. 

An  Envisu  R220  Spectral  Domain  OCT  (Bioptigen,  Durham,  NC)  system 
was  acquired,  along  with  a  lens  specific  to  the  focal  length  of  the  rat  eye.  An 
imaging  protocol  was  developed  to  generate  a  clear  image  of  left  and  right 
eyes,  each  centered  on  the  optic  nerve.  Rats  are  anesthetized  by  inhalation 
of  4%  isoflurane,  and  the  pupils  are  dilated  with  1%  tropicamide.  The  OCT 
images  (Fig.  8)  will  be  analyzed  using  either  in-house  thickness 
measurement  or  the  automated  Diver  2.0  software  from  Bioptigen. 


Figure  8.  Optical  coherence  tomography  (OCT)  is  used  to  monitor  changes  in  the  retinal 
thickness  (A).  A  series  of  retinal  images  are  taken  in  a  radial  pattern  around  the  optic  nerve 
of  the  eye  (B). 
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Visual  Evoked  Potential  (VEP) 

VEP  was  added  to  the  protocol  to  help  distinguish  between  eye  injury  and 
injuries  to  the  optic  nerve  or  visual  cortex.  Briefly,  eyes  are  exposed  to  a 
light  stimulus  and  resulting  visual  cortex  activity  is  measured.  The  light 
stimulus  is  made  up  of  10  white  LEDs  with  0.064  watts  per  LED.  Each 
stimulus  consists  of  a  half  second  duration  light  pulse  from  the  array 
followed  by  no  less  than  a  half  second  period  during  which  the  array  is  off. 
The  stimulus  is  repeated  15  times  for  each  eye.  The  brain  signal  from  the 
visual  cortex  is  measured  using  Ag/CI  electrodes  (Grass  Technologies,  RI). 
The  signal  is  collected,  fed  through  an  amplifying  circuit  (built  in-house),  and 
then  recorded  using  LabView.  The  signal  is  filtered  post-hoc,  segmented, 
averaged,  and  analyzed  in  MATLAB.  A  bandpass  filter  (1-100  Hz)  and  notch 
filter  (60  Hz)  are  used  to  isolate  the  desired  brain  waves.  The  program 
isolates  and  averages  the  15  repeated  signals  for  each  eye  to  eliminate 
portions  of  the  signal  not  directly  related  to  the  VEP  (Fig.  9).  The  primary 
metrics  recorded  for  the  VEP  are  the  peak  positive  and  negative  amplitudes 
as  well  as  the  delay  between  stimulus  and  peak  response.  Our  initial  studies 
with  VEP  resulted  in  a  low  signal  amplitude.  This  is  likely  due  to  the  effect  of 
the  anesthesia  on  the  brain  signal.  For  future  studies,  we  have  switched 
from  ketamine  to  isoflurane  as  the  anesthetic. 


Figure  9.  Visual  cortex  signals  fom  VEP.  Fifteen  signals  are  averaged  to  determine  the  brain 
response  to  visual  stimulation. 

son/  5;  Complete  histology  for  ocular  injury  from  blast  exposure  (Year  1  - 
Year  3) 

Following  euthanization  at  the  appropriate  time  points,  animals  are  perfusion 
fixed  with  10%  formalin.  The  brains  are  extracted  and  stored  in  10% 
forrmalin.  The  eyes  are  extracted  and  stored  in  a  mixture  of  2% 
gluteraldehyde  and  2.5%  formaldehyde.  Both  tissues  will  be  processed  for 
axonal  injury  and  ocular  trauma. 


8 


Aim  3:  Identify  changes  in  vitreous  protein  expression  that  correlate 
with  visual  system  injury 

SOl/1/  1:  Collect  vitreous  samples  following  every  animal  experiment  in  Aim 
2.  (Year  1-Year  3) 

The  eyes  from  half  of  the  animals  in  each  group  will  be  used  to  evaluate 
cytokine  biomarkers  of  ocular  trauma.  Currently,  eyes  have  been  removed 
from  animals  in  the  first  group  and  are  being  stored  in  a  -70°C  freezer  until 
the  protocol  for  protein  analysis  is  validated. 

SOl/1/  2;  Assay  samples  for  NfH  (a  marker  for  retina!  deterioration),  VEGF,  IL- 
10,  MCP-1,  and  MIPS.  (Year  1 -Year  3) 

The  protocol  to  perform  the  assays  is  currently  in  development.  Equipment 
purchase  for  sample  analysis  is  underway. 

KEY  RESEARCH  ACCOMPLISHMENTS: 

•  Acquired  shock  tube 

•  Characterized  pressure  profile  throughout  tube  both  experimentally 
and  computationally 

•  Added  animal  mount  to  shock  tube 

•  Determined  and  created  robust  filtering  protocols  for  pressure  signals 

•  Designed,  programmed  and  manufactured  vision  behavior  system  for 
rats 

•  Purchased  OCT  and  developed  a  successful  imaging  protocol 

•  Added  VEP  to  enhance  interpretation  of  results 

•  Build  VEP  system  and  developed  appropriate  filtering  protocol 

•  Began  blast  animal  studies 

•  Began  lOP  characterization  from  blast 

•  Received  animal  use  approval  from  University  of  Utah  and  Department 
of  Defense 

•  Received  human  study  IRB  approval  from  University  of  Utah,  Veterans 
Hospital  of  Salt  Lake  City,  and  Department  of  Defense. 

•  Developed  a  collaboration  with  the  polytrauma  group  at  the  VASLC  to 
obtain  access  to  medical  records  of  military  personnel  involved  in  a 
blast. 

•  Submitted  requests  for  medical  record  numbers  of  blast  victims 
treated  at  the  University  of  Utah  hospital. 

REPORTABLE  OUTCOMES: 

In  YEAR  1  of  project,  the  work  has  been  presented  at  two  conferences: 

•  International  Symposium  on  Computer  Methods  in  Biomechanics 
and  Biomedical  Engineering  (Appendix  A  &  B) 

•  Annual  Utah  Biomedical  Engineering  Conference  (Appendix  c) 
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The  work  was  also  presented  to  the  polytrauma  group  at  the  veterans 
hospital  in  Salt  Lake  City  (Appendix  D) 

CONCLUSION:  The  successful  completion  of  the  studies  proposed  in  this  4 
year  project  will  form  the  basis  for  understanding  the  temporal  and  chemical 
progression  of  visual  system  injury  following  blast  exposure.  In  the  first 
year,  all  the  infrastructure  and  product  development  was  completed  to 
successfully  achieve  the  stated  goals  of  the  study.  Clinical  and  experimental 
work  is  underway.  The  proposed  project  was  enhanced  by  custom-making 
the  vision  behavior  device  and  adding  visual  evoked  potential  to  the  battery 
of  tests.  Once  completed,  the  results  from  these  studies  will  lead  to  earlier 
diagnoses  of  visual  system  injury,  expand  our  understanding  of  the  time- 
dependent  response  of  the  visual  system  to  blast,  and  help  identify 
treatment  options  for  mitigating  vision  loss  from  blast  exposure. 

APPENDICES: 

•  Appendix  A  -  Abstract  accepted  to  the  International  Symposium 
on  Computer  Methods  in  Biomechanics  and  Biomedical  Engineering 

•  Appendix  B  -  Poster  presented  at  the  International  Symposium  on 
Computer  Methods  in  Biomechanics  and  Biomedical  Engineering 
Appendix  C  -  Abstract  accepted  to  the  9*^'^  Annual  Utah  Biomedical 
Engineering  Conference 

•  Appendix  D  -  PowerPoint  presentation  to  the  VA  Polytrauma  group. 
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APPENDIX  A  -  Abstract  accepted  to  11th  International  CMBBE 


Proceedings  of  the  International  Symposium, 
Computer  Methods  in  Biomechanics  and  Biomedical  Engineering 
April  3  -  7,  2013,  Salt  Lake  City,  Utah,  USA 


VERIFICATION  OF  A  SHOCK  TUBE  FINITE  ELEMENT  MODEL  FOR  BLAST  OCULAR 

INJURY  SIMULATIONS 


Daniel  F.  Shedd  and  Brittany  Coats 


Department  of  Mechanical  Engineering 
University  of  Utah 
Salt  Lake  City,  Utah,  USA 


INTRODUCTION 

Blast  exposure  is  a  leading  eause  of  eye  injury  for  the  US  Army 
[1].  Typieally,  oeular  injury  oeeurs  from  explosive  shrapnel  and 
debris,  but  reeently  many  soldiers  have  developed  vision  defieits  6-12 
months  following  a  blast  exposure  without  any  signs  of  injury  [2].  To 
better  understand  the  meehanisms  behind  blast  indueed  eye  trauma  and 
subsequent  vision  loss,  we  propose  to  develop  a  finite  element  (FE) 
model  of  an  eye  exposed  to  pressure  waveforms  assoeiated  with  blast. 
As  a  first  step  toward  this  model,  we  sought  to  develop  an  aeeurate 
representation  of  the  shoek  wave  pressure  profiles  assoeiated  with  real 
world  blast  exposure.  Therefore,  the  primary  objeetive  of  this  study 
was  to  develop  and  validate  the  shoek  tube  FE  model  to  be  used  in  the 
eye  blast  simulations.  A  seeondary  objeetive  of  this  study  was  to 
identify  the  shoek  tube  design  eriteria  required  to  generate  the  desired 
pressure  profiles  for  experimental  oeular  blast  studies  in  rats.  These 
experimental  data  will  be  used  to  validate  the  future  eye  blast 
simulations. 

METHODS 

To  simulate  the  pressure  profiles  in  the  shoek  tube,  an  Eulerian 
analysis  was  seleeted.  Helium,  hydrogen,  and  air  were  used  as  driver 
seetion  gases.  Air  was  used  in  the  driven  seetion  (Figure  1).  Eaeh  gas 
was  represented  using  the  ideal  gas  equation  of  state.  Material 
properties  used  in  the  model  are  listed  in  Table  1.  Initial  gas 
assignment  for  the  Eulerian  mesh  was  performed  using  referenee 
eylinders.  Temperature  was  used  to  set  the  initial  pressures  of  the  gas 
for  eaeh  of  the  tube  seetions.  The  open-ended  shoek  tube  was 
represented  as  a  rigid  body  with  an  outer  diameter  of  1.52  em  and  a 
length  of  4  m.  This  length  was  broken  up  into  a  1  meter  driver  seetion 
and  a  3  meter  driven  seetion.  The  effeet  of  gravity  was  eonsidered 
negligible.  The  membrane  separating  the  tube  seetions  was  not 


ineluded  in  the  model  beeause  the  time  period  investigated  was 


Figure  1.  Geometrie  and  Eulerian  seetions  of  the  finite  element 
model. 

The  Eulerian  spaee  was  meshed  with  linear  hexahedral  elements 
and  analyzed  in  ABAQUS/Explieit  (Dassault  Systemes,  Paris).  A 
eonvergenee  study  was  performed  using  a  driver  gage  pressure  of 
256.1  kPa  to  mateh  previous  simulations  [3].  Convergenee  of  peak 
overpressure  was  assessed  using  six  mesh  densities  (from  2500  to 
925,344  elements)  at  2.5  and  5  ms  post-rupture.  The  analysis  was 
performed  along  the  eentral  axis  of  the  driven  seetion.  Convergenee 
simulation  results  were  eompared  to  theoretieal  shoek  tube 
ealeulations  [4].  A  deviation  less  than  5%  from  the  theoretieal 
ealeulations  was  eonsidered  aeeeptable. 


Density 

I^rrt3 

OynamJc  Viscosity 
k/m-s 

Gas  Constant 

J/kE-K 

Specific  Heat 
J/kg-K 

Air 

1.171 

l.a5E-05 

2S7 

716 

Hefium 

0.16674 

1.9  6  E-OS 

2077 

3120 

Hydrogen 

0.0899 

S.7SE-06 

4126 

10135 

Table  1.  Ideal  Gas  properties  for  air,  helium,  and  hydrogen. 


Pressure  values  (integration  point,  gage  pressure)  generated  by  a 
200  kPa  initial  driver  overpressure  were  reeorded  along  the  entire 
length  of  the  eentral  axis  of  the  driven  seetion  from  0  to  12  ms  post¬ 
rupture  (Figure  3).  The  peak  overpressure  every  10  em  was  extraeted. 

Simulations  of  three  driver  gases  (air,  helium,  and  hydrogen) 
were  run  (Figure  4).  The  peak  overpressure,  peak  overpressure 
duration  (defined  as  pressure  >.95  peak  value),  and  time  eonstant 
(pressure  >.632  peak)  for  eaeh  driver  gas  was  reeorded  at  the  eentral 
axis  of  the  tube  at  the  longitudinal  midpoint  of  the  driven  seetion 
(1.5m  from  muzzle).  The  loeation  was  ehosen  to  allow  enough  spaee 
from  the  driver  seetion  for  the  shoek  wave  to  fully  form  while 
allowing  enough  distanee  from  the  tube  mouth  to  eolleet  meaningful 
data.  For  eaeh  driver  gas,  initial  temperature  eonditions  were  seleeted 
sueh  that  initial  driver  overpressure  was  200  kPa. 

RESULTS 

The  model  eonverged  to  within  1.5%  from  the  analytieal 
predietion.  The  density  of  the  final  Eulerian  mesh  was  338676  nodes 
and  312500  elements.  An  example  simulation  result  at  t=1.8  ms  is 
shown  in  Figure  (2).  Peak  overpressure  of  the  wave  deereased  as  the 
shoek  wave  traveled  down  the  tube  (Figure  3).  The  rate  of  the  deerease 
was  -4.7  kPa/m.  Hydrogen  gas  resulted  in  the  highest  peak 
overpressure  at  131.0  kPa  (helium:  124.4  kPa;  air:  92.0  kPa.)  and  the 
shortest  peak  overpressure  at  75.2  ps  (helium:  199.4;  air:  774.5).  The 
time  eonstants  of  the  pressure  waves  resulted  in  a  similar  trend: 
Hydrogen  .550  ms.  Helium  .687  ms.  Air  2.926  ms  (Figure  4). 


Figure  2.  Cross-seetion  view  of  overpressure  at  shoek  front,  t=1.8 
ms.  Data  was  reeorded  along  midline  elements. 

DISCUSSION 

Peak  shoekwave  overpressure  diminished  as  the  shoekwave 
traveled  towards  the  mouth  of  the  tube.  The  trend  was  -4.7  kPa/m, 
whieh  equates  to  a  loss  of  14.1  kPa  over  the  length  of  the  driven 
seetion,  or  10.4%  of  initial  driven  seetion  overpressure.  Shoek  tube 
theory  indieates  that  the  shoek  wave  overpressure  should  remain 
eonstant  throughout  the  driver  seetion  after  the  initial  30-50%  of  the 
tube.  Real  shoek  tubes  do  exhibit  diminishing  overpressure  along  the 
tube  due  to  boundary  effeets  and  energy  dissipation. 

Changing  the  driver  gas  impaeted  the  peak  overpressure,  peak 
duration,  and  time  eonstant  of  the  resulting  shoek  waves.  Light  gases 
(Helium,  Hydrogen)  behaved  eloser  to  ideal  gases  and  better 
approximated  the  Friedlander  waveform. 

The  environment  outside  of  the  shoek  tube  was  not  modeled. 
Therefore,  when  the  shoek  wave  nears  the  end  of  the  tube  the  aeeuraey 
of  the  data  eannot  be  ensured.  By  observation  of  trends  in  the  data,  the 
last  10%  (30  em)  of  the  tube  appears  to  have  been  affeeted.  This 
suggests  that  test  subjeets  used  to  develop  validation  data  for  the  future 
eye  blast  simulations  must  be  positioned  prior  to  this  loeation. 

Membrane  rupture  was  assumed  to  oeeur  immediately  prior  to 
simulation.  Real  membranes  rupture  by  splitting  into  flaps  of  material 
whieh  partially  impede  gas  flow  and  introduee  asymmetry  into  the 
shoek  front.  Membrane  seraps  ean  also  be  earried  down  the  tube  with 


the  shoek  front  or  the  subsequent  air  jet.  In  this  simulation,  rupture 
was  modeled  as  if  the  membrane  disappeared  upon  rupture. 

Future  studies  will  plaee  an  eye  model  within  the  simulated  shoek 
tube  and  traek  the  pressure  inside  and  outside  of  the  eye  over  time,  as 
well  as  the  strain  that  develops  within  the  eye.  These  data  will  be 
assoeiated  with  experimental  visual  dysfunetion  experiments  from 
blast  exposure.  Charaeterization  of  the  deerease  in  overpressure  along 
the  length  of  the  tube  will  be  used  to  determine  the  ideal  plaeement  of 
test  subjeets.  The  driver  gas  results  will  be  used  to  fine  tune  the 
magnitude  and  duration  of  the  blast  shoek  wave. 


Figure  3.  Pressure  time  profiles.  Curves  represent  10  em  of  driver 
length.  Note  the  deerease  in  peak  overpressure. 
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Figure  4.  Shoek  wave  at  midpoint  of  driver  seetion  for  air, 
helium,  and  hydrogen  driver  gases. 
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Introduction 


[Results 


Blast  exposure  is  a  leading  cause  of  eye  injury  for  the  US  Army 
[1].  Typically,  ocular  injury  occurs  from  explosive  shrapnel  and  debris, 
but  recently  many  soldiers  have  developed  vision  deficits  6-12  months 
following  a  blast  exposure  without  any  signs  of  injury  [2].  To  better 
understand  the  mechanisms  behind  biast  induced  eye  trauma  and 
subsequent  vision  ioss,  we  propose  to  deveiop  a  finite  eiement  (FE) 
modei  of  an  eye  exposed  to  pressure  waveforms  associated  with  biast. 
As  a  first  step  toward  this  modei,  we  sought  to  deveiop  an  accurate 
representation  of  the  shock  wave  pressure  profiies  associated  with 
reai  worid  biast  exposure.  Therefore,  the  primary  objective  of  this 
study  was  to  deveiop  and  vaiidate  the  shock  tube  FE  modei  to  be  used 
in  the  eye  biast  simuiations.  A  secondary  objective  of  this  study  was 
to  identify  the  shock  tube  design  criteria  required  to  generate  the 
desired  pressure  profiies  for  experimentai  ocuiar  biast  studies  in  rats. 


Methods 


To  simuiate  the  pressure  profiies  in  the  shock  tube,  an  Euierian 
anaiysis  was  seiected.  Heiium,  hydrogen,  and  air  were  used  as  driver 
section  gases.  Air  was  used  in  the  driven  section  (Fig.  1).  Each  gas 
was  represented  using  the  ideai  gas  equation  of  state.  The  open- 
ended  shock  tube  was  represented  as  a  rigid  body  with  an  outer 
diameter  of  15.2  cm  and  a  iength  of  4  m.  The  effect  of  gravity  was 
considered  negiigibie.  The  membrane  separating  the  tube  sections 
was  not  inciuded  in  the  modei  because  the  time  period  investigated 
was  immediateiy  subsequent  to  membrane  rupture. 

The  Euierian  space  was  meshed  with  iinear  hexahedrai  eiements 
and  anaiyzed  in  ABAQUS/Expiicit  (Dassauit  Systemes,  Paris).  A 
convergence  study  was  performed  using  a  driver  gage  pressure  of 
256.1  kPa  to  match  previous  simuiations  [3].  Convergence  resuits 
were  compared  to  theoreticai  shock  tube  caicuiations  [4].  A  deviation 
iess  than  5%  from  the  theoreticai  caicuiations  was  considered 
acceptabie. 

Pressure  vaiues  (integration  point,  gage  pressure)  generated  by 
a  200  kPa  initiai  driver  overpressure  were  recorded  aiong  the  entire 
iength  of  the  centrai  axis  of  the  driven  section  from  0  to  12  ms  post¬ 
rupture.  The  peak  overpressure  every  10  cm  was  extracted. 
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Figure  2.  Pressure  time  profiles.  Each 
curve  represents  history  output  of  single 
element.  One  curve  generate  for  each  10 
cm  along  the  driven  section.  Peak 
overpressure  decreased  as  a  function  of 
distance. 
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Figure  3.  Shock  wave  at  midpoint  of 
driver  section  for  air,  helium,  and  hydrogen 
driver  gases.  Each  curve  represents  the 
overpressure  at  one  element  across  time. 
The  peak  overpressure  increased  with 
decreasing  driver  gas  density. 


Figure  4.  Comparison  of  result  from  FEM 
simulation  and  results  from  experimental  shock 
tube  with  200  kPa  driver  pressure. 

Experimental  blast  produced  low  overpressure 
due  to  membrane  inefficiency.  Blast  duration 
was  accurately  predicted. 


Figure  5.  Cross-section  view  of  overpressure  (Pa)  of  the  shock  front  at  t=1.8 
ms.  This  simulation  used  a  helium  driver  gas  with  initial  pressure  set  to  200  kPa 
(29  psi).  Data  was  recorded  along  central  axis  elements  for  analysis  to  avoid 
edge  effects  near  the  boundary  of  the  tube. 


The  modei  converged  to  within  1.5% 
from  the  anaiyticai  prediction.  The  density  of 
the  finai  Euierian  mesh  was  338676  nodes 
and  312500  eiements.  Peak  overpressure  of 
the  wave  decreased  as  the  shock  wave 
traveied  down  the  tube  (Fig.  2).  The  rate  of 
the  decrease  was  -4.7  kPa/m.  Hydrogen  gas 
resuited  in  the  highest  peak  overpressure  at 
131.0  kPa  (heiium:  124.4  kPa;  air:  92.0  kPa, 
Fig.  3)  and  the  shortest  peak  overpressure 
at  75.2  ps  (heiium:  199.4;  air:  774.5).  The 
time  constants  of  the  pressure  waves 
resuited  in  a  simiiar  trend:  Hydrogen  .550 
ms,  Heiium  .687  ms.  Air  2.926  ms.  The  200 
kPa  simuiation  was  vaiidated  against  shock 
tube  data  from  a  200  kPa  biast  with  air  as 
the  driver  gas  (Fig.  4).  A  cross-section  view 
of  pressure  within  the  driven  section  is 


Conclusion 

Euierian 

Space 


Figure  1.  (A)  Experimental  shock  tube.  6"  internal  diameter.  Blast 
triggered  via  pneumatically  activated  arrowhead  rupturing  BoPET 
membranes  of  varied  thickness.  Instrumented  with  1  MS/s  pressure 
sensors  (PCB  113B26)  along  the  length  of  the  driven  section.  (B)  Geometry 
of  finite  element  model. 


The  Euierian  shock  tube  simuiation  was  successfui  in  capturing  the  deveiopment  of  a  Friediander  wave  as  the 
shock  wave  traveied  down  the  tube.  Additionaiiy,  distinct  differences  in  peak  overpressure  and  pressure  wave 
duration  couid  be  seen.  These  simuiations  wiii  assist  in  the  design  of  future  experimentai  studies  investigating  vision 
ioss  from  biast  exposure.  FEM  simuiations  were  aiso  abie  to  capture  the  time  duration  of  actuai  experimentai  data,  but 
overestimated  peak  overpressure.  This  is  iikeiy  due  to  the  effect  of  membrane  rupture  on  the  peak  overpressure, 
which  was  not  incorporated  into  the  modei. 

Future  studies  wiii  piace  an  eye  modei  within  the  simuiated  shock  tube  and  track  the  intraocuiar  pressure  (lOP) 
over  time,  as  weii  as  the  strain  that  deveiops  within  the  eye.  These  data  wiii  be  associated  with  experimentai  visuai 
dysfunction  experiments  from  biast  exposure.  Characterization  of  the  decrease  in  overpressure  aiong  the  iength  of  the 
tube  wiii  be  used  to  determine  the  ideai  piacement  of  test  subjects.  The  driver  gas  resuits  wiii  aid  in  seiecting 
appropriate  driver  gas  for  desired  overpressure  magnitude  and  duration  parameters.  This  modei  provides  a  tooi  for 
experimentai  design  for  biast  studies  using  shock  tubes  and  for  simuiating  the  effects  of  biast  on  bioiogicai  structures. 


Driver 

Gas 

Dynamic 
Density  Viscosity 
(kg/m3)  (kg/m-s) 

Gas 

Constant 

(J/kg-K) 

Specific 

Heat 

(J/kg-K) 

Air 

1.177 

18.5*10-6 

287 

716 

Helium 

.16674 

1.96*10-6 

2077 

3120 

Hydrogen 

.0899 

8.75*10-6 

4126 

10185 

Tabie  1.  Ideal  Gas  properties  for  air,  helium,  and  hydrogen  driver  gases. 
Air  was  used  as  the  driven  section  gas  for  all  simulations. 
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Appendix  C  -  Abstract  accepted  to  the  9'^'^  Annual  Utah  Biomedical 
Engineering  Conference 


Blast  exposure  is  a  leading  cause  of  eye  injury  for  the  US  Army.  Typically,  ocular  injury  occurs  from  explosive 
shrapnel  and  debris,  but  recently  many  soldiers  have  developed  vision  deficits  6-12  months  following  a  blast 
exposure  without  any  signs  of  injury.  To  better  understand  the  mechanisms  behind  blast  induced  eye  trauma 
and  subsequent  vision  loss,  we  propose  to  develop  a  finite  element  (FE)  model  of  an  eye  exposed  to  pressure 
waveforms  associated  with  blast.  As  a  first  step  toward  this  model,  we  sought  to  develop  an  accurate 
representation  of  the  shockwave  pressure  profiles  associated  with  real  world  blast  exposure. 
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APPENDIX  D  -  PowerPoint  presentation  to  the  VA  polytrauma  group. 
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Progression  of  visual  dysfunction  following 
blast  exposure 


Characterize  ocular  injuries 
associated  w/  and  w/out  TBI 
Identify  timelines  of  injury 
diagnosis/awareness 
Determine  early  identifiers  of 
visual  injury 

Funded  by  USAMRAA  Vision  Research  Program 
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Determine  blast  mechanics 
associated  with  ocular  injury 
Identify  timelines  of  ocular 
injury 

Determine  early  identifiers  of 
visual  injury 
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Clinical  Approach 

Step  1  -  Retrospective  review 

•  Identify  people  with  a  history  of  blast  exposure  and  visual  exams 

•  Record  complaints  of  vision  issues  and/or  diagnosed  issues 

•  Record  presence/absence  of  visual  degradation  of  disorders 

•  Record  time  after  exposure  that  visual  disorder  diagnosed 

Step  2  -  Prospective  study 

•  Enroll  people  who  have  been  exposed  to  a  blast 

•  Record  all  information  known  about  blast  (date  occurred,  distance  from 
source,  cause  of  blast) 

•  Perform  ocular  exams  (OCT,  optokinetics,  etc)  every  6  months  for  up  to  2  yrs. 

•  Record  self-assessed  vision  questionnaires  every  6  months  for  up  to  2  yrs. 

Funded  by  USAMRAA  Vision  Research  Program 
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